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To understand if sequential star formation is taking place 
at the periphery of the Hn region Sh2-219. We present 
12 CO J — 2 — > 1 line observations of this region, obtained at 
the IRAM 30-m telescope (Pico Veleta, Spain). In the optical, 
Sh2-219 is spherically symmetric around its exciting star; fur- 
thermore it is surrounded along three quarters of its periph- 
ery by a ring of atomic hydrogen. This spherical symmetry 
breaks down at infrared and millimetre wavelengths. A molec- 
ular cloud of about 2000 M lies at the southwestern border of 
Sh2-219, in the Hi gap. Two molecular condensations, elon- 
gated along the ionization front, probably result from the in- 
teraction between the expanding H n region and the molecular 
cloud. In this region of interaction there lies a cluster contain- 
ing many highly reddened stars, as well as a massive star excit- 
ing an ultracompact H n region. More surprisingly, the brightest 



^vq parts of the molecular cloud form a 'chimney' , perpendicular to 



< 



the ionization front. This chimney is closed at its south-west ex- 
tremity by Ha walls, thus forming a cavity. The whole structure 
is 7.5 pc long. A luminous Ha emission-line star, lying at one 
end of the chimney near the ionization front, may be respon- 
(N ' sible for this structure. Confrontation of the observations with 
models of H n region evolution shows that Sh2-219 is probably 
| 10 5 yr old. The age and origin of the near-IR cluster observed 
■ on the border of Sh2-219 remain unknown. 
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Abstract. Massive-star formation triggered by the expansion of H n regions. 



1. Introduction 

We have previously suggested (Deharveng et al. 2003a) that 
Sh2-219 is a case of sequential star formation. Our main 
argument was the presence of a deeply embedded cluster at 
the periphery of the Hn region, just beyond the ionization 
front. However, the molecular environment of Sh2-219 was 
unknown. Here we present high-angular resolution molecular 
observations, as well as new near-infrared (near-IR) and radio 
continuum observations, followed by a discussion of the 
morphology of the whole region in the light of these new data. 

Sh2-219 is an optically-visible Hn region, of diameter 
4.4 pc, lying at a distance of 5.0 + 0.8 kpc (Deharveng et 
al. I2003a> . in the direction I = 159?355, b = T.592. It is a 
nearly perfect Stromgren sphere around a central 09.5V excit- 
ing star (see Sect. 5). The mass of the ionized gas is 39 M Q 
(Leahy 1997 ). Sh2-219 is surrounded by a thick ring of atomic 
hydrogen, except for a gap in the southwest. This ring has a 
FWHM of 2-3 pc, a mean density of 9 atoms crrT 3 , and a mass 
of 97 M (Roger & Leahy 1993 1. We have previously reported 
(Deharveng et al. 2003a) the detection in the near- and mid-IR 
of a cluster southwest of Sh2-219, at the border of the ionized 
region, and in the direction of the H i gap. This cluster contains 
highly-reddened objects (Ay up to 14 mag for objects detected 
in all three JHK bands). It also contains an Ha emission-line 
star (no. 139 in Deharveng et al. 2003a) presenting a near-IR 
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excess, and an ultracompact (UC) Hn region (Leahy 1997). 
The coordinates of these objects are given in Table 1. 

Figure 1 shows Sh2-219 as it appears in the optical and the 
near infrared. 



Table 1. Coordinates of the objects discussed in the text 



Object 




RA (2000) 


Dec (2000) 




h 


m s 


o 1 It 


Exciting star of Sh2-219 a 


4 


56 10.61 


+47 23 35.4 


Ho- emission star no. 139 a 


4 


56 03.76 


+47 22 58.2 


UC H ii region 15 


4 


56 02.1 


+47 23 07.0 



a Deharveng et al. (2003a) 
b This paper 



2. Observations 

2.1. CO observations 

In December 2002 we observed the emission of the molecu- 
lar gas associated with Sh2-219, in the 12 CO J — 2 — > 1 line 
(v=230 GHz), using the IRAM 30-m telescope (Pico Veleta, 
Spain). We mapped an area of 14' x 15' with the HERA nine- 
channel heterodyne array (Schuster et al. 2004). The beam size 
of the telescope is 1 1'.'7 at this frequency. The elements of the 
array are arranged in a 3 x 3 matrix, with a separation of 24" 
on the sky between adjacent elements. The data were acquired 
by drifting the telescope in right ascension using the standard 
'on-the-fly' technique. A reference position was taken 10' east 
of the Ha star. It was checked to be free of emission down to 
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Fig. 1. Sh2-219. Left: composite colour image of Sh2-219 in the optical. North is up and east is left. The size of the field is 
5!1 x 5.'2. Pink corresponds to the Ho- 6563A emission, and turquoise to the [SII] 67 17-6731 A emission, enhanced near the 
ionization front. The arrow points to the Ha emission-line star (no. 139 in Deharveng et al. 1200 3a I. Right: composite colour 
image of Sh2-219 in the optical and the near-IR. Blue corresponds to the Ha emission of the ionized gas, and orange to the K' 
emission of the stars (new K' observations discussed in Sect. 4); the cross shows the northeastern limits of the K' frame. 



0.15 K in antenna temperature (ltr, in a velocity interval of 
0.25 km s~') between -70 and -6 km s _I . The beam pattern on 
the sky was rotated by 18.5 degrees with respect to the right 
ascension axis by means of a K-mirror mounted between the 
Nasmyth focal plane and the cryostat of the heterodyne array. 
When drifting the telescope in right ascension, two adjacent 
rows are separated by 7", which results in a slightly under- 
sampled map in declination (the Nyquist sampling step is 5'.'5). 
Details about the HERA array and the K-mirror can be found 
at http://www.iram.fr/IRAMFR/PV/veleta.htm 

A digital autocorrelator with a spectral resolution of 78 kHz 
was used as a spectrometer; the resolution was degraded later to 
obtain a velocity resolution of 0.25 km s at 1.3 mm. The ob- 
serving conditions were typical for the time of the year, with 
typical system temperatures of 400 to 500 K. Pointing was 
checked every 90 minutes by scanning across nearby quasars; 
it was found to be stable to better than 3". 

Supplementary observations in the C ls O J — 2 — > 1 
and /=l-»0 transitions at 219.560319 GHz and 
109.782182 GHz, respectively, were carried out at two 
positions in the nebula, using the 'standard' heterodyne 
receivers at the IRAM 30-m telescope. The weather conditions 
were very good, with system temperatures of 145 K and 260 K 
at 3 mm and 1.3 mm respectively. 

Given the extended size of the CO-emitting region, the an- 
tenna temperature is a satisfactory approximation to the CO 
line brightness. On the other hand, for the C ls O emission, 
which is much more compact, the main-beam brightness tem- 
perature scale r m b is a reasonable approximation to the intrinsic 



line brightness. We adopt values of 0.70 and 0.50 for the main- 
beam efficiency at the frequencies of the C ls O J — 1 — > and 
J — 2 — > 1 transitions respectively. 

In all the maps presented hereafter, the coordinates are ex- 
pressed in arcsecond offsets with respect to the position of the 
Ha emission-line star (Table 1). 

2.2. Other molecular observations 

As part of a larger observational program, we observed Sh2- 
219 with the Very Large Array on 2005 November 16. We 
searched for emission from methanol at 44 GHz, water at 
22 GHz, and ammonia (in the (1,1) and (2,2) lines) at 23 GHz. 
The array was in its most compact configuration, providing 
spatial resolutions of l'/5 in the 7 mm band and 3'.'3 in the 
1 .3 cm band. The pointing center was RA (2000)= 4 h 56 m 03 s .3, 
Dec (2000)= +47° 22' 57", and the primary beam was 1' and 
2', respectively in the 7 mm and 1.3 cm band. Further observa- 
tional details will be found in Deharveng et al. (2007, in prepa- 
ration). With 4<x detection limits of 45, 23, and 8 mJy beam -1 
for the methanol, water, and ammonia lines, respectively, we 
failed to detect any of these molecular transitions. 

2.3. Near-IR observations 

Near-IR observations were carried out through H and K' broad- 
band filters, and [Few] and H2 narrow-band filters, with the 
NICS camera (Baffa et al. l2"0TTfl at the 3.58-m Telescopio 
Nazionale Galileo (TNG, Canary Islands) on January 4, 2004. 
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The plate scale is 0725/pixel, yielding a field of view of 
4:2 x 4:2. The H 2 and [Fen] filters are centred at 2.12 /mi 
and 1 .67 //m respectively. Four pairs of on-source, off-source 
images on a dithered pattern were acquired through the H2 
and [fen] filters (total on-source integration time 10 minutes). 
Two series of four dithered on-source and off-source images 
were acquired through the K' filter (total on-source integration 
time 80 sec) and one series of on-source and off-source images 
through the H filter (total on-source integration time 80 sec). 
Flat-field images were obtained from sky observations at sun- 
set through all four filters. 

Each frame was first corrected for crosstalk using the rou- 
tine provided on the TNG web page (http://www.tng.iac.esi. 
Data reduction was then performed by using standard IRAF 
routines. After fiat-fielding all frames, the four off-source im- 
ages of each series were median-combined and subtracted from 
each on-source image of the same series in order to remove the 
background emission. The sky-subtracted images thus obtained 
were corrected for distortion and bad pixels, and were averaged 
together. 

Aperture photometry was done using the DAOPHOT 
(Stetson 1987) package in IRAF. We selected an aperture ~ 1 
FWHM in radius and the sky intensity was estimated from 
its median value in a ~ 1-FWHM wide annulus, between 2 
FWHM and 3 FWHM from the aperture centre (FWHM ~ 078 
in H and ~ 077 in K') . Stars were found using DAOFIND and 
the list was then corrected for false or missed detections after a 
visual check of the images. 

Although the seeing was good, the night was barely pho- 
tometric, so we calibrated our H and K' photometry by using 
that performed on the same field by Deharveng et al. ( 2003a 1. 
We calculated the difference in magnitudes for each coinciding 
source and plotted it against our derived instrumental H — K' . 
A colour term was present in ^'because of the difference of the 
K and K' filters. After a best-fit calibration, an rms residual of 
~ 0.2 mag was present in the difference between our photom- 
etry and that of Deharveng et al. ( 2003a 1, probably due to the 
different seeing, to the fact that we are comparing aperture and 
PSF-fitting photometry, and also possibly to the intrinsic vari- 
ability of young stars. Our resulting photometry is in the same 
HK system as Deharveng et al. (2003a). Sources brighter than 
~ 1 1 mag in H or K are in the non-linear range of the detector; 
in these cases, we adopted the measurements by Deharveng et 
al. ( I2003a> . The limiting detection magnitudes are ~ 18 in K 
and ~ 18.5 in H, much deeper than obtained by Deharveng et 
al. d2003at . We estimate a completeness limit of K ~ 16.5 by 
plotting the number of detected sources in 1-mag bins. A total 
of 616 near-IR sources with detection both in the H and in the 
K band and 53 with detection in only one band were found over 
the field. 

As for the narrow -band images, we removed the continuum 
contribution falling in the band by scaling both the narrow- 
band and the corresponding broad-band frame (K' for Hi and 
H for [Fe n] according to the signals measured for a few stars 
in the field (which are line-free, pure-continuum sources) and 
subtracting one from the other. The resulting frame contains 
only line emission. 



2.4. Radio continuum observations 

A radio map of Sh2-219 with a resolution of 40", obtained at 
20 cm with the VLA D array by Fich (1993), showed a sin- 
gle radio-continuum source. A later 20 cm map, obtained by 
Leahy (119971 ) with the VLA C array, had better resolution and 
showed the presence of an unresolved source south-west of 
Sh2-219. But no details were given in Leahy's paper about this 
ultracompact component. 

We used archival VLA data in order to assess the radio 
properties of both Sh2-219 and the UC Hn region to the south- 
west. The observations were made at 20 cm and 6 cm in the D 
and C configurations (programme AL216) and at 20 cm in the 
B configuration (programme AF346). The B array 20 cm data 
were consistent with the C array 20 cm data, and for simplicity 
will not be discussed further. 

The AL216 observations were made in January 1991 (C 
array) and April 1991 (D array). For both observation dates the 
primary calibrator was 3C48 and the secondary calibrator was 
0435+487. The adopted flux densities for 3C48 are 15.41 Jy at 
20 cm (1.46 GHz) and 5.16 Jy at 6 cm (4.86 GHz). The boot- 
strapped flux densities for 0435+487 were 1.25 and 0.52 Jy (at 
1.46 and 4.86 GHz, respectively) in January, and 1.24 and 0.51 
Jy in April. Data reduction followed standard procedures for 
VLA continuum data. 

The data obtained in the C configuration at 20 cm provide 
a uv-coverage similar to those obtained in the D configuration 
at 6 cm; hence these are the appropriate datasets for comparing 
the flux densities at the two frequencies, for both Sh2-219 and 
the UC Hn region. The C configuration's 6 cm data provide 
the highest angular resolution and hence give a stronger limit 
on the size of the UC H 11 region. The angular resolutions and 
noise levels are 14"xl2" and 140//Jy beam -1 (C configuration, 
20 cm), 15" x 1 1" and 65 //Jy beam -1 (D configuration, 6 cm), 
and 4" x 3" and 60 //Jy beam -1 (C configuration, 6 cm). The D 
configuration 6 cm image is shown in Fig. 2. 

In addition, a 1.3 cm (23.71 GHz) continuum image was 
produced from the ammonia observations mentioned in Sect. 
2.2. This image has a resolution of 276 x 273 and a noise level 
of 0.2 mJy beam -1 . 

3. The molecular gas 

3. 1 . Distribution and kinematics 

The distribution of integrated CO emission as a function of ve- 
locity is shown in Fig. 3. CO emission is detected to the south- 
west of Sh2-219 as a circular cloud of diameter ~ 12 pc, adja- 
cent to the H 11 region; the CO lines have typical antenna tem- 
peratures of 6-10 K. There is no molecular emission observed 
in the direction of the H 1 ring, nor is any CO associated with a 
possible photodissociation region surrounding the ionized gas. 

The overall CO emission covers a velocity range of ~ 
6 km s -1 . The molecular emission associated with Sh2-219 is 
particularly bright in the range -26.5 km s -1 to -25.0 km s -1 , 
where the bulk of the emission is detected; the systemic veloc- 
ity of the cloud is ~ -26kms -1 . This is in good agreement 
with previous velocity measurements (V(CO) = -24.5 km s -1 , 
Blitz et al. fT982l V(CO) = -25.0 kms -1 , Wouterloot & 
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Fig. 3. Channel maps of the CO J = 2 — > 1 emission as observed with HERA and integrated in velocity intervals of 0.5 km s ; 
the central velocity is marked in the upper left corner of each panel. The contour levels are 1, 2, 3, 5, . . . , 13 Kkms" 1 . The 
positions of the exciting star of Sh2-219 and of the Ha emission-line star (at 0,0) are marked with white crosses. 



Brand 1989). No emission is detected at -31.0 km s -1 , the ve- 
locity measured for the ionized gas by Fich et al. ( 1990). Note 
however that a mean velocity of -24.6 km s was measured by 
Georgelin & Georgelin ( 1970) for Sh2-219, in good agreement 
with that of the adjacent molecular cloud. 

The CO lines are relatively narrow, about 1 .5 km s _1 . There 
is no evidence for a bipolar outflow which could be the signa- 
ture of ongoing star formation in the cloud. However, detection 
of compact outflows, typical of UC H n regions, may be ham- 
pered by the spatial undersampling of the data and the low spa- 
tial resolution (one telescope beam corresponds to 0.3 pc at a 
distance of 5 kpc). 

Figure 4 shows the overall distribution of the ionized gas, 
of the atomic hydrogen, and of the molecular gas. These com- 
ponents are clearly separated. Two maxima of CO emission are 
found at the ionization front, at offset positions (-19", +9") 
and (+23", -32"). Figure 5 shows that these condensations are 
elongated along the ionization front. The former peak coin- 
cides with the UC Hn region (see Sect. 5); this condensation 
has a beam-deconvolved FWHM of 17" (0.4 pc). The latter 
peak coincides with some highly-reddened stars in the near-IR 



cluster. Its beam-deconvolved FWHM is similar (0.5 pc). Also, 
the CO emission map suggests the presence of a 'chimney' in 
the molecular cloud which nearly splits the cloud in two. This 
chimney is oriented southwest to northeast (approximate pro- 
jection angle of -135°). It is best detected in the blueshifted 
and redshifted gas (see e.g. the panels at -27 and -24kms~ 1 
in Fig. 3), but is nonetheless present in the ambient gas at 
-26kms~ 1 . The Ha emission-line star lies near the head of 
the chimney, and close to its axis. This chimney has an approx- 
imate length of 5 pc and a width of 0.5-1.0 pc. 

To gain more insight into the kinematics of the cloud, we 
have studied the profile of the CO emission along a southeast- 
to-northwest cut across the chimney, i.e. perpendicular to its 
main axis. Figure 6 shows the CO line intensity along the cut, 
as a function of the velocity. The gas appears to possess a 
small velocity gradient of 0.17 km s~' pc~' from Y = -100" 
to Y = +350". The cavity in the gas is visible as the region 
of low emissivity near Y = 0". Two components, with veloci- 
ties shifted with respect to that of the ambient gas, are detected 
near this position, one on the northern side at -24.0 km s _1 and 
one on the southern side at -27.5 kms -1 (Fig. 6). These com- 
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Fig. 5. Left: Composite colour image of Sh2-219. The Ha emission of the ionized gas appears in blue, the stellar K' emission in 
orange. The red contours are for the velocity-integrated temperature distribution of the CO emission. Note how the CO maxima 
are elongated along the ionization front. The arrow points to the Ha emission-line star. Right: The CO contours are superimposed 
to the radio continuum emission map at 20 cm. 
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Fig. 2. Radio continuum image obtained with the VLA at 6 cm. 
The angular resolution of the image (indicated in the lower left 
corner) is 15" x 11". The contour levels are multiples of the 
65 yt/Jy beam -1 image rms. 



ponents form the walls of the cavity. The southern wall reaches 
an extreme velocity of -31 kms -1 , and the northern wall of 
-22 km s . This is quite puzzling: it may indicate that the in- 
ner walls of the chimney are in rotation around the chimney 
axis, with the northern wall receding from us and the southern 
wall approaching us. 
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Fig. 4. The ionized, atomic and molecular components of the 
Sh2-219 complex. The Hn region (DSS-2 red survey) appears 
as a grey scale. The black contours show the ring of atomic ma- 
terial surrounding the ionized gas (from Roger & Leahy 1993 1. 
The thin red contours are for the velocity-integrated CO emis- 
sion. The first contour is 5.8 K km s _1 and the contour interval 
is2.9Kkms 1 . 



Faint CO emission at ~ 6 kms~' is observed south of 
Sh2-219 (A6 < 250"). This is most probably foreground emis- 
sion, unrelated to this region. 
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Fig. 6. Left: Distribution of the CO emission integrated between -28.0 km s _1 and -22.5 km s _1 ; contours are 0.5, 1, 2, 3 . . . times 
5.75 K km s _1 . The position of the Ha star at (0", 0") is marked with a white cross. Coordinates are arcsecond offsets with respect 
to the Ha star. The positions of the two C I8 observations are indicated by black dots. Right: Velocity-position diagram along a 
cut perpendicular to the chimney's main axis. The cut was made 28" southwest of the Ha emission star, along the south-east to 
north-west direction shown as a dashed line in the left panel. Y=0 corresponds to the projection of the Ha star on the cut. 



3.2. Physical conditions 

The brightest CO emission is detected at the head of the chim- 
ney, in the two condensations adjacent to the ionization front, 
with antenna temperatures of up to 12 K, indicating a gas ki- 
netic temperature of ~ 15 K. The CO antenna temperature 
along the chimney is lower, typically 8-10 K. 

The molecular material associated with Sh2-219 has been 
mapped only in 12 CO (2-1), so for the mass determination 
we rely on the empirical relation between cloud mass and in- 
tegrated emission j T(CO)dv. We use the CO to Ho conver- 
sion factor X - 1.9 x 10 20 mol cm~ 2 (K km s -1 ) -1 , as de- 
termined by Strong & Mattox (1996 1. Including all emission 
above 1 K km s" 1 (~ 3.5cr) this yields M C o = 2100 M Q . The 
equivalent radius of the cloud above that emission level is r e = 
y/Aiea/n = 5.8 pc. The linewidth of the composite cloud spec- 
trum is A V — 2.4 km s , but this reflects the fact that different 
parts of the cloud have different velocities (cf. Fig. 3), rather 
than turbulence. Using the typical profile width (1.5 km s _1 ; see 
Sect. 3.1) instead, and assuming a density profile oc r~ 2 , we find 
a virial mass M vir = 126 x r e (A V) 2 = 1800 M , in good agree- 
ment with the CO mass derived above, and consistent with the 
cloud as a whole being in gravitational equilibrium. 

Likewise, the mass of the gas associated with the chimney, 
corresponding to j T( n CO)dv > 10 K km s"\ is 770 M Q ; 
the equivalent radius of this emission region is about 2 pc (only 
an indicative value of course, as the chimney emission region 
is elongated). The composite CO spectrum of the chimney 
emission has an equivalent width of ~ 3.1 km s _1 - wider than 



that of the cloud as a whole, because the chimney gas shows 
the largest range in velocity (see Figs 3 and 6). Here also, the 
individual profiles have a width of about 1.5 km s . 

Two positions in the cloud, indicated in Fig. 6, were ob- 
served in the C 18 J — 2 — > 1 and J — 1 — > lines. The lines 
are narrow (1 km s -1 ) and weak: typically T m \, = 0.2 K. At 
offset position (-70", -10") - behind the northern border of 
the chimney - we find T^ h = 0.36 K and 7"™ = 0.22 K. 
At this position, the CO J — 2 — > 1 antenna temperature is 
« 7 K. We have estimated the physical conditions from a 
simple calculation in the large-velocity gradient approxima- 
tion, assuming uniform density and temperature in the emit- 
ting region. Adopting a kinetic temperature of 15 K and a 
linewidth of 1.0 km s _1 , we obtain n(H%) = 8.0 x 10 3 cirT 3 and 
N(C 18 0) = 2 x 10 14 cirT 2 . Adopting a standard C 18 abun- 
dance of 2 x 10~ 7 (Frerking et al. 1982 1, this corresponds to 
an H 2 gas column density N(H 2 ) of 1.0 x 10 21 cirT 2 . The C ls O 
lines are optically thin with t 21 = 0.054 and t 10 = 0.015. The 
excitation temperature of the J — 2 — > 1 transition is only 1 1 K, 
i.e. it is not thermalized. On the other hand, the J — 1 — > tran- 
sition is close to thermalization. We also observed the C ls O 
emission at the offset position (+30", -84"), near the southern 
border of the cloud. The lines intensities are comparable and 
we derived a similar density and a similar column density in 
the gas traced by C ls O. 

If we assume that the emission detected at the two observed 
positions is representative of the bulk of the chimney gas, the 
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derived NQlz), together with the 2 pc radius found from the 
12 CO-observations, leads to a mass of 270 M Q for the chimney 
gas. This (very uncertain) value is almost three times less than 
what is derived from the 12 CO-emission. If the two CO clumps 
are density enhancements (as opposed to regions with an el- 
evated temperature), they would contain much more material 
than the walls of the chimney and would consequently exhibit 
stronger C ls O lines, which would lead to a higher mass esti- 
mate. As no C 18 observations were carried out towards the 
two clumps this cannot be confirmed. 

4. The stellar population of the cluster 

A colour-magnitude (K vs. H-K) diagram for all sources de- 
tected in both bands is shown in Fig. Most of the stars fall 
in a strip within Ay ~ 10 mag from the ZAMS, but a few 
sources exhibit larger extinctions (and, possibly, near-IR ex- 
cesses). This confirms the results of Deharveng et al. (2003ai: 
their colour-colour diagram (J - H vs. H - K; see their fig. 4) 
shows that most of the sources are consistent with reddened 
main sequence stars (with Ay < 15 mag), whereas only a 
few exhibit a near-IR excess. However, our present data un- 
veil the presence of a larger number of near-IR sources with 
H-K > 1 mag than are seen in Deharveng et al. (2003a). These 
may be the youngest members of the cluster. Note that the 
brightest of these reddened sources is star no. 1 39 of Deharveng 
et al. (2003a) - the Ha emission-line star - which is ~0.6 mag 
brighter in K than given by their photometry. This may indicate 
that this source is variable. 

To better examine the cluster morphology, we have mapped 
the stellar surface density by counting all sources (above the 
completeness limit) in a 20" x 20" window with sampling ev- 
ery 10" both east- west and north-south. The result is shown 
in Fig. 8a. The locations of the sources are marked either by 
small filled squares (H - K < 1) or large squares (H - K > 1). 
The surface density distribution appears slightly elongated in 
a direction roughly parallel to the ionization front of the H 11 
region. The reddened sources (presumably also the youngest 
stars in the cluster) are even more clearly aligned along it. Most 
of these sources are located towards the two most prominent 
CO clumps, at offsets of (-19", +9") and (+23", -32"). 

The radial distribution of the stellar surface density can 
be shown by counting all sources (above the completeness 
limit) in 10"-wide annuli centred at the cluster's centre, at 
RA (2000)= 4 h 56 m 03*9, Dec (2000)= +47° 22' 51". The 
result is shown in Fig. 8b. All derived parameters for the clus- 
ter are listed in Table 2. From the radial density distribution 
of Fig. 8b, we have estimated the number of cluster mem- 
bers (down to the completeness limit) as in Testi et al. ( 1997). 
The cluster size is derived from the radial density distribution; 
Table 2 gives the radius at which the surface density falls al- 
most to the background value (~15 stars arcmkr 2 ). The com- 
pleteness limit in mass, M comp i, can be estimated from the com- 
pleteness limit at K by adopting a set of evolutionary tracks. 
From Palla & Stahler ( 1999), K = 16.5 corresponds to a pre- 
main-sequence star 10 6 yr old of 0.4 Mq or 10 7 yr old of 
1.1 Mq. Assuming an Ay of 10 mag, these mass limits increase 
to 1 .0 Mq and 1 .4 Mq, respectively. Hence the cluster popula- 
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Fig. 7. K vs. H — K for all sources with detection in both bands. 
The dashed lines mark the completeness limit, whereas the thin 
line is the locus of the ZAMS, for Ay = and a distance of 
5 kpc; the absolute magnitudes are from Allen ( 1976 ), and the 
colours from Koornneef d!983> . The thick lines are isochrones 
for pre-main-sequence stars 1 Myr in age (with Ay = and 10 
mag, left to right) with masses from 0. 1 to 6 Mq (from Palla & 
Stahler 1999 ). The length of the dotted reddening line originat- 
ing from a B0 star corresponds to a visual extinction of 30 mag; 
the adopted extinction law is from Rieke & Lebofsky ( 1985). 
The positions of the exciting star of Sh2-219 (09.5V) and of 
the He emission star (no. 139) are given. 



tion is probably well sampled down to stars of ~ 1 Mq. This 
means that the number of members down to K = 16.5 is a 
lower limit to the total number of members; assuming a Scalo- 
like IMF (Scalo 1998), the correction amounts to a factor be- 
tween two (for M comp i = 0.4 M Q ) and six (for M comp i = 1 Mq), 
without including unresolved companions and stars below the 
hydrogen burning limit. 
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Fig. 8. a Stellar surface density distribution (stars arcmin~ 2 ); coordinates are arcsec offsets from the Ha emission-line star no. 139. 
Small squares mark the positions of sources with H - K < 1, and large squares mark the positions of sources with H - K > 1. 
The locations of the UC Hn region and of star no. 139 are indicated by crosses, b Radial stellar density distribution, centred at 
the surface density peak. 



Table 2. Cluster parameters 



Number of members with K < 16.5 ~ 175 

Fraction of members with H - K > 1 22% 
Radius 1.5 pc 

Expected number of members 3 down to 0.1 Mq 380 (1040) 
Expected number of members 3 with M > 10 Mq 2 (4) 
Total stellar mass" 280 (750) Mq 

a assuming all the stars down to 0.4 (1) Mq have been detected 



If we knew the cluster age we could compare it with the 
dynamical age of the Hn region. As shown in Fig.0 a sim- 
ple overlay of pre-main-sequence isochrones indicates that the 
stellar population is compatible with pre-main-sequence stars 
10 6 yr old. However, the cluster age cannot be constrained 
by the magnitude-colour diagram, as any isochrones between 
zero and several megayears are compatible with it. Information 
about the age comes from the number of sources showing a 
near-IR excess. Deharveng et al. (2003a) found only a few 
sources with a near-IR excess in the J - H vs. H - K diagram. 
Even assuming that all our sources with H - K > 1 do exhibit 
a near-IR excess, this amounts to only 22% of cluster members 
(see Table 2). This shortage of sources with a near-IR excess 
may then suggest a cluster where star formation started more 
than 10 6 yr ago (Haisch et al. l2001> . 

As noted above, the distribution of sources with H - K > 1 
suggests that most recent star formation occurs parallel to 
the ionization front, possibly in the two most prominent CO 
clumps. This is confirmed by the fact that the north-western 
clump hosts a UC Hn region and that the Ha emission-line 
star lies at its edge. As for the south-eastern clump, our Hi im- 



age (Fig. 9) shows faint elongated line emission in its direction. 
This could be the signature of a jet from a young star; we have 
tentatively identified the driving source with the star indicated 
in Fig. 9 (K = 13.7, H - K = 1.3). But this could be simply 
fluorescence emission from the surface of the molecular cloud, 
illuminated by UV radiation leaking from the H n region. No 
other knots of H2 line emission have been found in the field 
and no [Fe 11] line emission has been detected. 

5. The ultracompact Hh region 

Figure 2 gives the 6 cm radio continuum map of Sh2-219, 
showing the presence of an ultracompact Hn region at its 
southwestern border. Leahy ( 1997 > showed the 20 cm coun- 
terpart of this map, i.e. the 20 cm C configuration map of VLA 
program AL216, while we show the 6 cm D configuration map. 
The focus of Leahy ( 1997 1 was the Sh2-219 region; that paper 
did not discuss the UC H 11 region to the south-west. 

Sh2-219 has a flux density of 140 ± 7 mJy at 20 cm and of 
127 ± 6 mJy at 6 cm. These values are consistent with an opti- 
cally thin H 11 region, with an ionizing photon flux of log iVLyc = 
47.5. According to the calibration of Martins et al. (2005 1, this 
flux corresponds to an 09.5V exciting star, in good agreement 
with the K magnitude of the central exciting star of Sh2-219 
(Fig. 7). 

The UC Hn region has a flux density of 3.2 + 0.3 mJy at 
20 cm, 6.1 + 0.4 mJy at 6 cm, and 4.2 ± 0.8 mJy at 1.3 cm. 
To within the uncertainties, these flux densities are consis- 
tent with a thermal spectrum. The deconvolved source size is 
< 076 x 075, corresponding to < 15 mpc x 12 mpc. The min- 
imum ionizing photon flux required to maintain ionization of 
the region is log 7VL yc > 46.1, corresponding to a B0.5V or 
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Fig. 9. Enlargement of the #2 image of the cluster: left, line plus continuum emission within the band, right, after subtraction of 
the continuum emission. The arrows mark the faint emission feature, tentatively identified as a jet (right panel) and the possible 
jet-driving source (left panel). The Ha emission-line star lies at (0,0). 



B1.5V exciting star (according to Panagia fl 9731 and Smith et 
al. 2002 respectively). Using the 6 cm flux density and the up- 
per limit to the size, we derive a lower limit for the rms elec- 
tron density in the UC HII region of 1.3 X 10 4 cm 3 . This elec- 
tron density is slightly higher than the H2 density found in the 
molecular clumps; this is to be expected as the CO data are 
averaged over a much larger volume due to the lower spatial 
resolution. 

The density and temperature of the molecular clump sur- 
rounding the UC H 11 region are too low to provide a pressure 
high enough to confine the ionized gas. So it seems very un- 
likely that this UC H n region is old and has reached pressure 
equilibrium with its surroundings. Furthermore, the radius of 
the initial Stromgren sphere of a star with log A^yc = 46.1 in 
a medium of density 1.3 X 10 4 cirT 3 is about 13 mpc - almost 
exactly the upper limit we report for the size of the observed 
UC H 11 region. Thus it appears that the UC H 11 region is quite 
young, and has not yet had time to expand from its initial size. 

6. Discussion 

6.1. The Ha emission-line star no. 139 

In the following, we try to estimate the luminosity of the Ha 
emission-line star no. 139, observed at the southwestern border 
of Sh2-219; this star is identified in Fig. 1. This is the brightest 
K star of the near-IR cluster. 

The integrated Ha flux of Sh2-219, from spectrophotomet- 
ric observations (Caplan et al. 2000), is 2.5 x 10~ n erg s^cirT 2 . 
The Ha flux of the emission-line star was determined as fol- 
lows. First, all stars were removed from the Ha image using 
DAOPHOT, and the integrated emission of the H 11 region was 
then measured, exempt from stellar pollution, in the instrumen- 
tal units. This allowed us to calibrate the instrumental units 
(ADUxpixels) with respect to flux units (erg s^'crrT 2 ). Next, 
with DAOPHOT, the stellar fluxes were measured in instru- 
mental units from both the (original) Ha image and the [S 11] 



image. For all except the emission-line star the measured sig- 
nals in these two filters were proportional (since they mea- 
sure continuum emission at neighbouring wavelengths but with 
different filter widths). For the emission-line star, the mea- 
sured Ha minus the Ha expected from the [S 11] filter mea- 
surement gave us the star's instrumental emission-line sig- 
nal corrected for its stellar continuum. Using our calibration, 
this difference was then converted to an Ha emission flux of 
2.9 x 10~ 14 erg s^'cirT 2 . At a distance of 5 kpc this corresponds 
to an Ha luminosity of 8.6 x 10 31 erg s _1 (2.3 x 10~ 2 L Q ). 

The Ha emission-line star presents a strong near-IR excess 
which prevents us from accurately determining - and correct- 
ing for - its extinction. Assuming that its J - H colour is not 
affected by the near-IR excess, and that its intrinsic / - H is 
0.0 mag, typical of early type stars, the observed J - H of 
1.6 mag indicates a visual extinction of ~ 15 mag. This is close 
to the maximum visual extinction obtained for nearby stars 
in the cluster from measurements in the three ( JHK) bands 
(Deharveng et al. 2003a). The assumption that star no. 139 
has a visual extinction of 15 mag leads to an Ha luminosity 
of 6 x 10 36 erg s _I (~ 1600 L Q ); this is an upper limit, as its 
J-H colour is probably affected by the near-IR excess. A lower 
limit is obtained assuming a visual extinction of 4 mag, the 
extinction of interstellar origin affecting Sh2-219 (Deharveng 
at al. 2003ai; hence an Ha luminosity of 1.7 x 10 33 erg s _1 
(~ 0.45 Lq). The Ha luminosity of star 139 is typical of pre- 
main sequence objects with mass loss and - possibly - outflows 
(LevreaultQ][EEK 

The Ha star no. 139 is probably a high-mass Herbig Ae/Be 
star, similar to MCW1080 or V645Cyg (which display mass 
loss and are associated with outflows; Levreault 1988 1. It is 
also similar to object 25 in the AFGL4029 cluster (Deharveng 
et al. 119971 ). an emission-line object with an Ha luminosity of 
3.4 x 10~ 2 L o (uncorrected for extinction), strong extinction and 
a near-IR excess. 
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6.2. The CO chimney - Ha cavity 

Fig. 10 is a deep Ha image of Sh2-219, resulting from the addi- 
tion of six frames (total exposure time of one hour) obtained at 
the 120-cm telescope of the Observatoire de Haute-Provence. 
Most of the stars have been subtracted using DAOPHOT; the 
intensity of the filaments has been enhanced using Photoshop. 
This Ha image shows very thin low-brightness filaments south 
of Sh2-219. Some of these filaments are oriented east-west. 
They most probably belong to the foreground supernova rem- 
nant G160.9+2.6, also known as HB9 (Van den Bergh et 
al.[R73), whose distance is ~ 1.2 kpc (Leahy & Roger 1991 1. 
Other filaments, indicated by arrows, are oriented north-east to 
south-west and seem to form the walls of a cavity, closing the 
molecular chimney at its south-western extremity. For compar- 
ison, Fig. 10 also shows the molecular emission superimposed 
on the Ha image. The molecular chimney and the Ha cavity 
have the same axis. The length of the cavity is 7.5 pc, from the 
'head' of the chimney in the north-east, near star no. 139, to 
its southwest extremity. Large-scale Herbig-Haro jets or flows, 
driven by low-luminosity sources (Reipurth et al. U997l and ref- 
erences therein), are known to exist. The largest of these is 
HH 111, which has a total extent of 7.7 pc, and is driven by 
a 25 L Q class I source. Also, a few large-scale jets are driven 
by luminous sources. HH 80-81, for example, with a total ex- 
tent of 5.3 pc, is driven by a 2 X 10 4 L Q young B star (Marti et 
al. Qj393). However, none of these jets or flows are associated 
with molecular chimneys or cavities. 

The Ha emission-line star no. 139, located at the head of 
the chimney and near the axis of the cavity, is well positioned 
to be at the origin of such a structure. Near-IR spectra of this 
object are needed to ascertain its nature. 

6.3. General discussion 

Sh2-219 is an almost perfectly spherical Hn region around a 
central exciting star. Hence this region has probably evolved 
during most of its lifetime in a rather spherically-symmetric en- 
vironment. The ionized gas is surrounded by an H i photodis- 
sociation region (PDR), except in the south-west. In Sh2-219 
the density of the ionized gas is in the range 55-170 cm" 3 
(Deharveng et al. 2003a I, and, as in all H n regions, the electron 
temperature is of the order of 10 4 K. The surrounding atomic 
gas is of lower density, 9 cm" 3 , and of lower temperature, 100- 
500 K (Hosokawa & Inutsuka 2005 ). Thus the pressure in the 
ionized zone is higher than in the neutral surrounding zone; 
the Hn region is presently still in expansion, the ionization 
front eating away the atomic ring from the inside. According 
to Roger & Leahy dl993> the H i ring is not expanding; thus 
it is probably in pressure equilibrium with a lower density and 
partially ionized external medium. 

Using the very simple model of Hn region evolution by 




Dyson & Williams ( 1997 ), we calculate that if Sh2-219 formed 
and evolved in a uniform medium of density 10 2 cm" 3 , its 
present radius of 2.2 pc indicates a dynamical age of 1.5 x 
10 5 yr. This age estimate is very uncertain, because the as- 
sumption of evolution in a strictly uniform medium is unre- 
alistic. On the one hand, the exciting star was probably formed 



Fig. 10. The CO chimney - Ha cavity south-west of Sh2-219. 
The contours of velocity-integrated CO emission are super- 
imposed on the Ha image of Sh2-219. The first contour is 
5.8 K km s" 1 and the contour interval is 2.9 K km s" 1 . To high- 
light the Ha filaments, most of the low-brightness stars have 
been subtracted from the original frame, and the filaments were 
digitally enhanced. The arrows point to filaments forming the 
walls of the Ha cavity closing the molecular chimney. 



in a much denser molecular core; during that early period the 
region would have evolved more slowly than suggested by the 
simple model. But on the other hand, we know that the H n re- 
gion is expanding into a lower density medium and is therefore 
presently evolving more rapidly. 

Sh2-219 is surrounded by a large atomic photodissociated 
region. Roger & Dewdney flp92| and Diaz-Miller et al. ( 1998 1 
have presented models of the time-dependent evolution of pho- 
toionized and photodissociated regions formed by stars of var- 
ious effective temperature embedded in molecular material of 
uniform density. These models do not include hydrodynamics. 
They show that the photodissociated region is relatively more 
extended with respect to the ionized region when the exciting 
star has an effective temperature below 30000 K (i.e. for B0- 
B4 stars), and that it can persist a long time in a low-density 
medium. This applies well to Sh2-219. However these models 
are static (they assume that a dissociation equilibrium has been 
reached); the H2 reformation time is very long in a low den- 
sity medium, much longer than the dynamical timescale. Thus, 
as Diaz-Miller et al. (1998 1 have strongly emphasized, these 
models are not adapted to regions evolving in a low density 
medium, as is the case of Sh2-219. We therefore cannot apply 
these models to Sh2-219 to estimate its age from the relative 
sizes and masses of the ionized and atomic regions. 



Hosokawa & Inutsuka (2005 1 present new models where 
the dynamical evolution of an Hn region and of its PDR is 
analysed by calculating the radiation transfer and the thermal 
and chemical processes with a time -dependent hydrodynamics 
code. They illustrate their models with the case of an 06 
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star evolving in a 10 3 cm 3 uniform medium, which applies 
well to Sh2-104 (Deharveng et al. l2003b> but not to Sh2-219. 
Hosokawa et al. (in preparation) have improved these models 
by considering evolution in a non-uniform medium, and by 
taking into account photodissociation by the far-UV Galactic 
background radiation. Hosokawa et al. have considered the 
case of Sh2-219. Their model assumes that its exciting star 
formed in a cloud composed of a small high-density core 
(0.07 pc, 10 5 cm~ 3 ) surrounded by a region of radially de- 
creasing (oc r~' 5 ) density. The model accounts for the size of 
the ionized region and of the photodissociation region, and 
for the densities in the ionized and atomic material. The size 
of Sh2-219 implies an age ~ 10 5 yr. The CO molecule is 
completely photodissociated in the PDR; this explains why 
no CO emission is detected around Sh2-219 (except in the 
south-west). 

We can imagine two widely different scenarios to account 
for the overall morphology of this region and its star formation 
history. 

- Scenario 1: The exciting star of Sh2-219 formed inside the 
near-IR cluster, and consequently has the same age. It was 
subsequently ejected from the cluster following dynami- 
cal interactions (for example this may result from a three 
body encounter implying a binary system and a third star; 
Kroupa^OOO 2001 1. If the age of the cluster is 10 5 yr 
(10 6 yr), the present location of the exciting star with re- 
spect to the cluster implies an ejection velocity > 20 km s _1 
(> 2 km s _1 ). Several points argue against this scenario: i) 
why is it the most massive star of the cluster which has 
been ejected? ii) why is the present radius of the Stromgren 
sphere (which depends on several independent parameters) 
equal to the distance covered by the star since its ejection? 
iii) the star is presently surrounded by a spherical H n region 
and by an almost spherical atomic PDR; how is this spher- 
ical symmetry compatible with an ejection (a highly non- 
isotropic process)? iv) the observed density of the ionized 
gas is higher than the observed density of the surrounding 
atomic material; if the star has been ejected from its dense 
parental molecular cloud into a low-density medium, how 
can we account for this? Also, it is difficult to eject a mas- 
sive object from a cluster. It requires, first, mass segrega- 
tion, which leads to a high density of massive objects in the 
very centre of the cluster. Ejection of massive objects gen- 
erally happens several megayears after cluster formation, as 
shown by Kroupa (2000 his fig. 6). 

- Scenario 2: The exciting star of Sh2-219 formed first, in 
a molecular condensation close to - but distinct from - 
the cloud presently observed. It formed an Hn region 
and photodissociated most of the surrounding molecular 
material. During its expansion the Hn region interacted 
with a nearby pre-existing molecular cloud, forming a 
layer of shocked compressed material at the surface of the 
cloud; this material forms the two molecular condensations 
presently observed at the southwestern border of the Sh2- 
219 H n region. Hosowaka et al.'s model predicts that, if the 
average density of the molecular cloud is 8 x 10 3 cirT 3 , the 



time scale of the layer fragmentation is comparable with the 
estimated age of Sh2-219. Therefore triggered star forma- 
tion seems possible, the observed cluster being a second- 
generation cluster whose formation was triggered by the 
expansion of the Sh2-219 Hn region. It contains one mas- 
sive star exciting the UC H n region. A Herbig Be star - the 
Ha emission-line star, or another still-hidden object - sub- 
sequently modified the structure of the adjacent molecular 
cloud, its strong wind carving a chimney through it. Here 
again several points argue against this scenario: i) we are 
unable to prove that the cluster is younger than the exciting 
star of Sh2-219. Indeed, the absence of maser emission and 
the small fraction of stars presenting a near-IR excess seem 
to indicate that the cluster is possibly older than 10 5 years; 
ii) here again the spherical symmetry of the ionized region 
is difficult to reconcile with the interaction, in one direction, 
with a nearby molecular cloud. 

What observations would allow us to choose between these 
two scenarios? A measurement of the proper motion of the ex- 
citing star of Sh2-219, to confirm or deny the ejection, would 
be most helpful; but this measurement is not presently possible, 
this region lying too far away. Deeper near- and mid-IR frames 
of the cluster would be helpful, allowing an estimate of its age 
via the fraction of objects with near-IR excesses, and via the 
presence or absence of deeply embedded (class or class I) 
objects. Spitzer IRAC observations of this region would prob- 
ably give this information. High resolution observations of the 
two CO clumps adjacent to the Hn region could confirm and 
give details about the interaction between the expanding H n 
region and the molecular cloud. 

7. Conclusions 

The molecular CO observations presented here have shown the 
presence of a molecular cloud adjacent to the Hn region and 
filling a 'gap' in the ring of atomic hydrogen surrounding the 
ionized gas. This cloud has a mass of about 2000 M Q . A most 
remarkable aspect of its morphology is a chimney, ~5 pc long, 
crossing the entire cloud. The Sh2-219 Hn region lies at one 
extremity of this chimney. The other extremity is closed off by 
filamentary Ho- walls. We suggest that a bright Ha emission 
star, located near the ionization front, may be at the origin of 
the strong wind or flow which has carved out this chimney. 

The Sh2-219 Hn region is in interaction with this adjacent 
molecular cloud. Two molecular condensations are observed at 
the interface between the H n region and the cloud. A near-IR 
cluster is observed in the direction of this interaction region. 
It contains highly reddened stars, a luminous Ha emission star 
which is also a near-IR excess object, and at least one early B 
star exciting an unresolved ultracompact H n region. However 
no CO outflow and no maser emission are observed in the di- 
rection of the cluster, suggesting that it is not very young. 

Sh2-219 is one of the Hn regions we proposed as can- 
didates for the collect and collapse process of star formation 
(Deharveng et al. 2005 1. We now know that star formation does 
not result from this process, as the observed molecular cloud 
probably pre-existed, and was not formed of material collected 
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during the expansion of the H n region. No completely satisfy- 
ing explanation has been found to explain massive-star forma- 
tion in this region. 
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